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ABSTRACT 


feapoineipnalemoede Of failure of Structural components in 
mechanical systems is fatigue. One method of predicting the 
Meobability of fatigue failure of a structural component is 
momcetermine the probability that the calculated cumulative 
fatigue damage index is greater than the critical damage 
index at failure. The cumulative fatigue damage index is 
represented as a random variable, and the critical damage 
mex 1S represented by the statistical variance of existing 
experimental data. A FORTRAN computer code using this failure 
@ereeriad is presented, which calculates the probability of 
meure £Or a structural component in the high-cycle fatigue 
regime under a random stress response environment, using both 
the Weibull and log-normal statistical distribution models. 
The Weibull model has been found to be the more conservative 
Medel in the low probability of failure region, which is con- 
Sistent with failure predictions between the two models using 


Meemelassical failure criteria of cyclic life. 
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Pee UD NeeeouUCT TON 


Fatigue 1S a many facetted phenomenon, influenced by such 
Gmvyerse factors as component surface condition, mean and maxi- 
Mees tress, prestrain, temperature, loading rates, stress con- 
Semerations, corrosion, axial orientation of loading, loading 
sequence and random stresses or strains incurred [Ref. 1] 
through [Ref. 5]. All of these factors, as well as the mate- 
meme roperties of the component undergoing the fatigue process, 
have an enormous statistical variation in their determined mag- 
mrmeddes. This varlation in parametric values has led to the 
meeelopment of reliability analysis techniques for predicting 
m@enuseful life of a component. 

Perlyeprediction models were developed from extensive case 
history data files. From this information, mean stress and 
cycles to failure were readily available, so the most preva- 
lent analysis models are based on cyclic stress history tech- 
meeues. 16 apply these models to a new design, extensive 
testing must be accomplished to obtain sufficient data to get 
mea-onaple resolution of predicted results. The testing prob- 
m@ebecomes even more Severe for components operating in the 
“infinite life" region of stress-life (S-N) curves. 

Many structural components (i.e., heat exchangers, off- 
Shore structures, etc.) are subject to random loadings, such 


as flow excitation, with generally small vibration amplitudes 


al 





G@aestress/strain cyclic lives in excess of ioe 


cles marr Or 
mis situation, the fatigue problem is one of high-cycle ran- 
dom fatigue [Ref. 6]. Due to a lack of experimental data and 
mee histories for structures in this regime, fatigue failure 
prediction models based on the cumulative damage index of the 
Palmgren-Miner Linear Damage Law are being proposed. 

The purpose of this research was to develop a FORTRAN com- 
puter code based on the cumulative damage index that could be 
Meeadeto predict the probability of failure for a structural 
component in the high-cycle fatigue regime under a random 
stress response environment. The computer code developed 
calculates the probability of failure for such a component, 
Memlizing both a Weibull and a log-normal statistical distri- 
bution procedure. The computer code was tested by using ex- 
[rine experimental data, with the predicted results of the 
two distributions compared to each other. Additionally, the 
trends of the predicted values are compared with published 
descriptions of the behavior of these two distribution models 
as used in cyclic life analysis. 

The predicted probabilities of failure for the two models 
were compared for both single component structures, and simple, 
evciple component structures. in the multiple member case, 
the weakest link analogy was used as the structure failure 
@riteria. In this analogy, the structure is assumed to fail 


when the weakest component of the structure fails. 


eZ 





[irsomaecument COltalns the analytic development of equa- 
tions used in the computer code, as well as extensive graphical 
presentation of calculated results. The computer code and a 
users guide are presented as appendicies. Recommendations for 
future work, such as incorporation of multiaxial fatigue analy- 


sis and effects of maintenance procedures are also given. 





II. BACKGROUND 


A. FATIGUE 

Fatigue has been described as a process of progressive 
localized permanent structural changes in a material due to 
Remains at some point in that material, which may result in 
@eeexs or fracture after a sufficient number of cyclic fluc- 
mmeerons {Ref. 7]. Fatigue can be classified into two cate- 
gories; low-cycle fatigue and high-cycle fatigue. In the low- 
cycle fatigue regime, plastic strain predominates and ductility 
controls performance. In the high-cycle fatigue regime, elastic 
strain dominates and strength control performance. In the high- 
Cycle fatigue regime, most of the fatigue life 1s spent in the 
Seeck initiation phase, with a relatively short portion of life 
spent in the crack propagation phase. 

imueemdat the analytic technique developed here contains mate- 
rial fatigue strength as one of its control variables, and util- 
izes stress vice strain analysis, the model is only applicable 
to high-cycle fatigue life predictions. No applicability to 
low-cycle fatigue is intended or implied, thus the following 
discussion is limited to high-cycle fatigue. Many factors 
affect fatigue life in the high-cycle fatigue regime. Two of 
these factors are solely stress-dependent functions, and deserve 
Some amplification. Additionally, a brief discussion of S-N 


curves is presented. 
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if Loading. Sequence 


The loading sequence effect has received extensive 
investigation. One simple, but descriptive presentation of 
maeewerfect was tendered by Dowling [{Ref. 8]. Figure 2.1 is 
fmeeoresentation of the potential orders of stress loading 


sequence. 





meomGe 2.) SEQUENCE EFFECTS ON CRACK INITIATION (from Ref. 8) 


For the high-low stress loading sequence (Fig. 2.l{(a)), 


a crack may initiate in the high stress range, with resultant 


failure in the low stress range after fewer operating cycles 


LS 





than would normally be expected. For the low-high loading 
sequence (Fig. 2.1(b)), application of a large number of cycles 
meee 1Ow Stress level does not significantly affect the number 
of cycles required for failure at the high loading level. For 
multiple changes in loading sequence (Fig. 2.1(c)), a crack 
amnitiated at a high level can propagate in the low level. In 
all three cases, cracks are assumed to initiate in the high 
Stress loading level, but to propagate in either stress loading 
ieee tnese Cracks are only cracks due to loading sequence, 
feoeeco not include cracks that would have developed regardless 
of loading level (i.e., constant amplitude stress loaded struc- 
tures can develop cracks). Thus, loading sequence has a signi- 
ficant effect on fatigue life. 
2. Mean Stress 

Realistic fatigue life prediction in cases where the 
mean stresses are large, relative to the fluctuating stresses, 
must account for the affect of the mean stresses. Due to the 
Mespataively low fluctuating stress amplitudes in the high-cycle 
fatigue regime, mean stress effects are very important. If 
the strain levels are high enough to cause repeated plastic 
meen, mean stresses are rapidly relaxed and mean stress 
effects can be neglected [Ref. 6]. However, high-cycle fatigue 
1s dominated by elastic strain, not plastic strain, so mean 


stress effects must be accounted for. 
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Various investigators have proposed models which 
estimate the effect of mean stress, o, on fatigue life. Most 
of these models are based on the following relationship: 


SCN = ace. - (N) (eam (22). 1) 


q 


where Ao(N) stress range including the effect of mean 

stress as a function of cycles, N 

AD ag (N) = equivalent stress range that would cause 
failure in the absence of a mean stress 
expressed as a function of cycles, N 


a(o, ) = mean stress correction coefficient 


The four most frequently used models for a(o,) are 
fae cally depicted in Figure 2.2. Examination of these 
Smeyres demonstrate a Significant difference in predicted safe 
values for the design stress range. Comparison of the models 
with experimental data for ductile steels indicates the Goodman 
and Soderberg curves are too conservative, while the Gerber 
parabola seems to be the best representation. The Goodman 
model does correspond rather well with brittle steel data, 
milemas tracture strength approaches ultimate strength, the 
SAE model approaches the Goodman models predicted value 


fer. 6]. 


Ih 





STRESS RANGE, Ado 


SODERBERG 





MEAN STRESS, Oo 


Figure-2.2 PROPOSED MODELS FOR THE EFFECT OF MEAN STRESS 


where oT 


(amon Rete 9) 


Goodman model: is COT) 
Soderberg model: 1 - Gil 2p) 

: _ 1 
SAE model: 1 CONE ) 


, : 2 
Gerber model: 1 (or ee 


mean stress 

ultimate tensile strength 

yield strength 

patigue Strength coetficient (or true fracture 


strength) 
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Peeeotress-Lifte Curves 

Stress-life (S-N) curves have been used for determin- 
Meee fatigue analysis for decades. Typically, S represents 
the applied stress, and N represents number of cycles, or life 
to failure. These plots usually are presented on semilog or 
log-log coordinates. Most S-N curves have a continuously slop- 
mieetrace, but some curves (primarily low strength steels), 
have a discontinuity, with an essentially horizontal segment 


: Ee TO GhiGhes. 


mom 10 
This marked cessation of non-zero slope led to the 
first definition of endurance (fatigue) limit. The fatigue 
emt is the limiting value of the median fatigue strength as 
N becomes very large. The median value is used due to the 
Me@eter in test values. Commonly, test data at several stress 
levels is obtained, and then S-N curves are drawn through the 
Mmeaman points, and thus represent fifty percent expected fail- 
ures [Ref. 10]. This statistical scatter, and subsequent 
median value representation is an added reason for using prob- 
Bpalastic fatigue analysis, vice deterministic fatigue analysis. 
Pataeue 1S usually classified as either low cycle or 
high cycle fatigue, with the transition occurring around ius 
cycles. Most high cycle fatigue test data terminates around 
iad cycles, and the preponderence of S-N curves in the litera- 


0 7 


fUre range from 10 to 10° cycles. This study is concerned 


with fatigue lives well above the published curves, i.e., 


lie, 





fatigue lives greater than io = 


Soles somal, extrapolation 
method of available data and curves was needed. 

Recent literature classifies fatigue lives in the 
range of interest as being in the "ultra-high cycle fatigue" 
region, rather than the usual phraseology such as fatigue life 
in the "infinite life’ region of the S-N curve, because of the 
ews Or poorly defined endurance limit of most materials. 
Manson [Ref. 11] presents an excellent discussion on all three 
ranges of fatigue life, and of importance to this study, pre- 
Somes a summary on ultra-high cycle fatigue extrapolation. 
This summary is based on the work of G. R. Halford, of the 
ASME Nuclear Pressure Vessel Piping Code Activity. In that 
reference material is available in which techniques for, and 
feet cication for, extrapolating known S-N curves into the 
moron Of interest in this study, the assumption of the avail- 
ability of the appropriate S-N curves for fatigue life predic- 


mrons 1S valid. 


Eee ALMGREN-MINER HYPOTHESIS 

In 1945 Miner developed a cumulative linear fatigue damage 
model, based on criterion suggested by Palmgren in 1924 [Ref. 
Meaemevichin the context of this hypothesis, fatigue failure 
occurs when the cumulative fatigue damage index reaches one(l). 
Because of its simplicity, the model has been widely used in 
fatigue analysis. However, it does not include consideration 


Baethe statistical variability existing in fatigue behavior. 





ii@emecritical value (i.e., the actual value of the cumulative 
damage index at failure, not the assumed value of 1.0 at 
memlure) is not always close to one (1), but has a range of 
memes from 0.18 to 23, with only a small portion close to one 
feet. 12]. 

Momreccount £Or the statistical variability of the critical 
damage index (denoted as A), two probability distribution models 
have been used. These are the Weibull and the log-normal dis- 
mmrpietions. Wirsching [Refs. 12, 13] suggests using the log- 
Mounal distribution for A with a mean of 1.00 and a coefficient 
@eaariation of 0.65, based on the collection of experimental 
data of several investigators. Buch [Ref. 14] also examined 
numerous sets of existing experimental data, investigating 
m@@eevariation of the critical damage index A , for different 
loading spectra and structural components. 

Selecting thirty-seven sets of data from Buch's paper, 
Oder statistics [Refs. 15,16] were used to reduce the data 
mig plot it on probability paper. A best fit analysis was 
conducted, and the Weibull distribution with a mean of 0.9 
meea coefficient of variation of 0.67 was found to be the 
best representation of the data, as shown in Figure 2.3. The 
Weibull distribution shape parameters b and c were determined 
momoe 08.9791 and 1.3289 respectively. 

From the mean and coefficient of variation of A , the 
Standard deviation of A was determined to be 0.4781. Accepted 


MmuGbabtitstic theory indicates 99.7% of all possible values of 


Zee 





any distribution will be contained within the range of the 
mean plus and minus three standard deviations. Thus the 
Seatistical range of the A's from Buch's data is 0.00 to 
2.33, excluding the tail regions, while the actual range is 
Meieeto 5.40. This 1s a significant deviation from the value 
1.0 of the Palmgren-Miner hypothesis, and demonstrates the 
need for the development of a probabilistic fatigue life pre- 


Gmecton model. 


fee ROBABILITY MODEL OF A 

The thirty-seven sets of A and the resultant shape para- 
meters, mean, standard deviation and coefficient of variation 
were used to test the developed Weibull distribution model. 
Mimeesame mean and coefficient of variation were used in the 
log-normal distribution so that a legitament comparison could 
be made, with the understanding that the data fit the Weibull 
fjetribution better than it did the log-normal distribution, 
so some discrepancy was expected. The mean and coefficient 
of variation would be the same, regardless of which distribu- 
tion actually applied, but tail region scatter affects may be 
different for each distribution. 

ieee Werbull Distribution 


The probability density function BOS is defined as: 


ae 


Cc Cc 
f(x) =e [El  expl-(2 1 (eqn 2.2) 
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Tiewcumulative distribution function eS) is defined as: 


Cc 
Betx) = 1 - exp [-(§) | (eqn 2.3) 


tere xX > 0, b> 0, andc > Q. 
The relationship between the shape parameters b and c, and 
the mean (x) and the standard deviation (9,) of the random 


Variable x are: 


1 


Aoo= ia (i+ =| (eqn 2.4) 
and 
O Spe GS] Saya tiees (eqn 2.5) 
x © G j 
= a Beet. sacl 
where IT (a) = Gamma function = J e ie dit. 


From equations 2.4 and 2.5, the coefficient of variation, CV 


of r {1 + 4 
cv= ~ =f ——}— -1 (eqn 2.6) 
x 2 


ie 


- Log-normal IDs ere allo hes olal 


The other, more frequently used, distribution function 
in fatigue analysis is the log-normal distribution function. 
It is more prevalently used due to the availability of a closed 
form mathematical formulation. The probability density func- 


tion £ (x) is defined as: 





1 Meee se 
a = exp[- (5) (BLS) | ; (eqn 2.7) 


nN 


Xx 


where x is the mean of In(x) and Gi, is the standard deviation 
of In(x). The cumulative distribution function BIS is defined 


as; 


In(x)-x 


0 
Xx 


2 
F(x) = exp [-7] Cie x 20 (eqn 2.8) 


alr 


- © 


or 


PC) = & co aa (eqn 2.9) 


where @ {[ ] is the standard form of the normal distribution 
mametions. For a random variable x, if x follows a log-normal 
@eweribution, then In(x) follows a normal distribution [Ref. 17], 
mamseethe closed form normal distribution functions can be used. 
This alleviates the necessity of using numerical integration 
routines to solve the distribution functions. 

The mean, x, and the standard deviation, oj, of random 


Variable x are: 


oe Ne 


X = exp [x + = ] (eqn 2.10) 


hes 





and 


Solving for x and o, in terms of x and o,: 


and 


O 
X 


O 
Xx 


ras 
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i 2 2 
exp [2x + of ] { exp [oy ] 


= iis | = ~ ln ct Bo) 


Y-In(x*) + Hie = pas ) 
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III. FORMULATION OF ANALYSIS 


A. INTRODUCTION 

In formulating a probability based analysis of fatigue 
life for structural components, the measure of fatigue fail- 
ure can be expressed as the probability of failure, P,. The 
fagbure criterion to be used can be defined such that the 
design cumulative fatigue damage index (denoted as D), is 
greater than the critical damage index A , at failure. Both 
mamaee indexes are subject to statistical scatter, thus the 
probability of failure is a two variable function. The prob- 


Bomaty of failure can be defined as: 


oO 


a Ds) = | Fy (x) £p (x) dx (eqn 3.1) 
or 
Pe = ]1- P..[D<4] =a | FH(x)£, Oc) dx ceciio2) 
where Fy» PA -weorcnecunmilative Gistnribution functions of D 
and A respectively. 
fh» fA = the probability density functions of D and 


A respectively. 
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Due to the paucity of experimental data in the cyclic 
stress range of interest, a time history simulation model 
must be used to generate a representative data base for input 
to the probabilistic analysis model. The simulation model 
used in conjunction with this study was developed by Y. S. 
Pameiket. 6]. In that a substantial portion of the developed 
source code consists of this model, pertinent aspects of the 
model are presented in conjunction with the development of the 


probabilistic analysis model. 


B. CUMULATIVE FATIGUE DAMAGE 

As previously discussed, the most widely used criterion 
for determining fatigue failure for variable stress loads is 
the Palmgren-Miner hypothesis. The cumulative fatigue damage 
index, D, is the sum of the damage incurred at each stress 


Mewel, and can be expressed as: 


fa 


(eqn 5.5) 


js! 
Ze, 
> 


where i = elapsed number of cycles at stress level s. 


N, = total number of cycles at stress level s. 
inducing failure in the component 
D, = damage incurred at stress level s; 


For continuous systems, D can be expressed by: 


; Not, CS) 
D = == GSO ds (eqn 3.4) 
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or 


wh 


At 


ii 
Doe Nae Intsy! (eqn 3.5) 


ere N, = the total number of cycles in the design life 
£,(S)ds = fraction of stress peaks between stress 
level s and level stds 


N(S) 


number of cycles at stress level s 


expected value (probabilistic mean) of Y 


Ey | 


Poemre, D > A, and from equations 3.4 or 3.5, the total 


member of cycles to failure (Ne) can be calculated as: 


ON 


c. 


Ne a Oa (eqn Seb) 


N- = —>—— (eqn 3.7) 


RANDOM VIBRATION ANALYSIS 


"A stochastic process is employed to analyze the vibration 
responses statistically, and to determine the type of pres- 
sure (e.g., narrow/wide band process) and the peak envelope 
distribution with statistical parameters (e.g., mean, stand- 
ard deviation). High-cycle fatigue test data for materials 
under random vibration are scarce and the commonly avail- 
able sinusoidal fatigue S-N curves are resorted to as a 
basis to evaluate the random fatigue life." [Ref. 6] 


For a given fatigue life curve, equation 3.6 can be eval- 


uated if ce!) the probability density function (PDF) of a 
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peak stress envelope is known, or can be approximated. In most 
cases, the PDF cannot be easily derived. However, if the stress 
hastory follows a Gaussian process with zero mean stress, closed 
meen sSOlutions for US) are available. 


fee Probability DSS Te be eis ese selene Peak Stress Envelopes 


iitesmnregilarity factor, a , can be defined such that: 


ao 
a= — (eqn 3.8) 
m 
O 
where the, = expected rate of zero (mean stress) crossings 
il expected rate of peaks 


If a approaches zero, a very peaky process is indicated, and 
a Gaussian distribution for £ (8) can be assumed. Therefore, 


eal ) ean be expressed as: 


Z 





£.(S) = ——- exp [- 2 ] 3; -7°<S<o (eqn 3.9) 


where oe is the standard deviation of the stress amplitude. 
If q approahces one, a narrow band process is implied, and a 
Rayleigh distribution for eles can be assumed. Thus, Peas) 


can be expressed as: 


2 
S) 


20, 





mcs) = —, exp[- = AO <2 (eqn 3.10) 
S S 
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When the cyclic stress history follows the stationary 
Gaussian process with zero mean stress [Ref. 18], the PDF for 


the peak stress envelope can be expressed by: 





—— 2 
ms) = 1-ca’ exp [- S ] + 
2 JOT 26 = (ete) 2 ) 
TO. S 
we 1 + Sea a } exp[- ° ] 
20 P ae: (eqn 3.11) 
S Of 7 ae 
> @ 


where erf(x) is the Gaussian error function, and can be 


Beepressed as: 


erf(x) = e dt (eqn 3.12) 


ate 


Z. Number of Stress Cycles at Stress Level Sj to 
Cause Failure: Nj 





Mie SOUrec dara for this analysis can be obtained from 
existing sinusoidal S-N fatigue curves, constant amplitude test- 
ame, Or stress limits specified in codes, such as the ASME 
mmeremee Vesse] and Piping Codes. For example, if S-N data 


Mmemavallable such that 


b 
Se = Te N. (eames. 3) 

then 5 -B 
Ne isl (eqn 3.14) 





where 


2 = one (eam o..13) 


In which b is the fatigue strength exponent of the particular 
material being analyzed. These values have been experimentally 
determined and tabulated for selected engineering alloys 

[Ref. 19]. 


Substituting equation 35.14 into 3.7, 


B 8 
ep | PIS 


A = N ——— 
Cr LCP hs 


N 


ile L (eqn 3.16) 


For a given stress history the expected value 5S, Ellice can 
be estimated from: 
Bp ere elie 8 ; 
ES | = = - S5 fears <1 7 } 


Eee CLIC STRESS/STRAIN ANALYSIS 

When stress histories follow other than a stationary 
Gaussian process, or have a none-Zero mean stress value, 
eeced formesolutions are not avallable for the PDF. The 
use Of the random vibration approach to determine E[S?] LS 
eeeeremely difficult. To alleviate this problem, cycle count- 
ing algorithms are used. These algorithms are very useful 
if cyclic stress/strain histories for the component are known. 


The major disadvantage of using cycle counting methods 1s 





that a tremendously large amount of data is required to obtain 
reliable fatigue life predictions. To compensate for the lack 
Sraeddita, a typical stress/strain time history data block is 
estimated [Ref. 20]. The generated data block is assumed to 
mepeat itself ad-infinitum, so that fatigue life prediction 
calculations for high-cycle fatigue can be accomplished. 
me cyele Counting Algorithms 
The object of cycle counting methods is to compare 
meow), irregular load histories with S-N curves developed 
from uniformly repeated simple load cycles, which are readily 
available. 
'All good counting methods must count a cycle with 
the range from the highest peak to the lowest valley and 
seek to count other cycles in a manor that maximizes the 
Mamees that are counted. This rule can be justified 
either by assuming that damage is a function of the mag- 
nitude of the hysteresis loon, or by considering that in 
merrcue (as in many other fields) intermediate fluctua- 
tions are less important than the overall differences 
between high points and low points." 
Phi cood counting methods count every part of 
every overall range once and only once. They also count 
smaller ranges down to some predetermined threshold once 
and only once. Three counting methods that achieve this 
Gbjective are well documented in the literature: range- 
mais, Lainflow, and racetrack." ([Ref. 10] 
All three methods can be used with, or without mean stress con- 
Siderations. The racetrack method is considered the most suit- 
able for condensing actual load histories. The rainflow method 


is the most popular method, and can be used with load histories, 


or strain histories [Ref. 10]. 


5:9 





mee Raintlow Cycle Counting Method 


Diesramrlow cyele Counting method is used in this study, 
because it can identify the stress range associated with low 
frequency components, and the linear damage law summation pro- 
cedure allows addition of damage index values from large and 
small stress ranges [Ref. 6]. 

The ''rainflow counting'' method is so named because it 
reminded its developers, M. Matsuishi and T. Endo [Ref. 21], 
Of rain flowing down a series of pagoda roofs. Rules are 
imposed on this downward flow so that cycles and half cycles 
can be differentiated and counted. 

In application, a sample stress time history, X(t), 
memeonverted to a point process of peaks and troughs as shown 
in Figure 3.1, with the peaks assigned even numbers, and the 
meouens assigned odd numbers. The time axis is oriented ver- 
tically, with the downward direction being the positive direc- 
meon. A “rainfiow'' is started at each peak and each trough. 
When a rainflow started at a trough comes to a peak, the flow 
is terminated if the magnitude of the opposite trough is less 
meerne Originating trough (e.g., Fig. 3.1{b); paths 1-4, 
moon 7-8 and 9-10). For a path originating at a peak, flow 
is terminated by a peak of greater magnitude than the origin- 
ating peak (e.g., paths 2-3, 4-9). If the rain flowing down 
a path intercepts flow from a previously counted path, count- 


ing of the current path 1s terminated so that each path is 
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meee o. i) EXAMPLE OF RAITNFLOW CYCLE COUNTING METHOD 
(from Ref. 22) 
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omieeecounted once (e.g., paths 3-3a, 6-6a, 8-8a). A new path 
M@eeee cycle count) is not started until the path under con- 


Sideration is terminated by one of the above rules. 


Eee t BULL DISTRIBUTION MODEL 

The probability density function £ (x) and the cumulative 
aeseribution function BOS) were defined as equations 2.2 and 
eee the probability of failure, Pe, was defined as equations 
leon) 5.2. For this study, equation 3.1 was chosen. Sub- 


Seenting A and D into Ee 4 6X) and fap) respectively, the 


new cumulative distribution function can be expressed as: 
e 
F(x) = 1 - exp 5 a (eqn 3.18) 
mmdeetie new probability density function can be expressed as: 
Chl Ch 


f(x) = “D c= So) (eqn 3.19) 
me bp op NP FD a 


Weibull distribution shape parameters for the 


where Dy» Cy 


CEltleal damage Index 


Dp? €p Weibull distribution shape parameters for the 
design damage index 

From these equations, the probability of failure can now be 
expressed as: 


Pie = | PF, (x) E(x) dx (eqn 3.20) 


e) 
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Solution of this integral requires the use of numerical 
integration techniques. The program developed uses the NONIMSL 
subroutine DQSF, which is a combination of Simpson's rule and 
Newton's three-eighths rule. This routine is a fourth order 
mmecrsion algorithm with fifth order truncation error, and was 


deemed to be accurate enough for the purposes of this study. 


F. LOG-NORMAL DISTRIBUTION MODEL 
Beet and D follow log-normal distributions, Iin(A) and 
meeetollow normal distributions. Therefore, define random 


variable Z as: 


Bec) = In(aA)-- In(D) (eqn 3.21) 


Assuming that A and D are independent random variables, the 


mean and the variance of In(Z) are expressed by: 


Z=A- D Fean S222 
and 
= es OD (eqn 3.23) 
where Z, A , D = mean of 1n(Z), In(A), and 1n(D) 


OF» Of; mee tance or lm), haGweeand in D) 
ve Peeler GD) ; In(A)] = P_tin(Z) < 0] 


1 oF t2 a 
= —— exp[-+]dt =o [--*] (eqn 3.24) 
Jot | 2 Os 


— 





The log-normal distribution shape parameter is simpler to 
obtain than the Weibull shape parameters. For the log-normal 
model, it is simply the ratio of the means of the two random 
Meeraples. Within this study, the shape parameter will be 
called the central safety factor. The central safety factor 


(n) is expressed as: 


n = (eqns +25) 


ot |e 


Deame equations 2.12 and 2.13, (A,0,) and (D, ip) can be 
expressed in terms of the central safety factor and the coef- 


Mmeerents of variation CV, amie CV py aS 


/ KR? 











A = Lin (ae . ino + A*) = se ecmimer2 ©) 
SC 

gt = ¥ - In(B?) + In(o,?+ D2) = ¥ In[1+CV,7] — (eqn 3.27) 

A _ = ve NH 2 

D = 1n(D?) - 5 In(c,? = B%) = in / 2 (eqn 3.28) 
1+CV)* 


on = -1n(D*) = In(op* + Dy) ae) In[1+CV)*] (eqn 3.29) 


S@imerttuting equations 3.26 through 3.29 into equations 3.22 


pid 5.23, 
= 1+CV,° 
een) + in (eqn 3.30) 
LON 
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and 
ee) in (i + CV,*) (1 + CV,°)] (eqn 3.31) 


Menee, equation 3.24 can be written as ([Ref. 25]: 


icv, 
Wiel) aes tea ag x 


ee 
Yin[(1+CV,*) (1+CV,)] 


Ceome.5..02)) 


Where ® [ ] is the standard argument for normal distribu- 
tions, as previously discussed. The program developed uses 
the IMSL subroutine MDNOR, which evaluates the absolute value 
Of the standard argument and returns the probability of sur- 
Vival. The probability of failure is then defined as one 


minus the probability of survival. 


G. WEAKEST LINK ANALOGY 

In the weakest link analogy, a multiple member structure 
momassumed to fail when its weakest member fails. That 1s, 
the remaining members cannot take on, without failing them- 
selves, the load born by the failed member. For a multiple 
component structure, the probability of failure can be 


defined as: 


Pe eee) 0 21 = 1,2,3,...,N] (eqn 5.5.5) 


Ree, 





oO 


a ~ 
me PC [i $I (kg SMBs 3 Ueto? oe eee N| 6 (eqn 63.54) 


Where A; = the critical damage index of the 1-th member 


the design damage index of the i-th member 
iy [(A; - D;) < 0 ] 1s the occurrence of at least 
one event (A, - D,) < 0 for any "i" 
N = the number of members 
iti [(A, 


ce 


- ey > 0 |] is the occurrence of all events 
> Dell &  ssiere ahi ah 
Assuming "4," and "D," are mutually independent, and the design 


damage index is the same for all members (Dd: ee nen 


oO 


N 
ee = | - Il ewe rer tx) 1 f(x) dx (eqn 3.35) 
i aa De D 
1=1 al 
where FA: (x) SGuMurative aGistribution funetion of A; 
al 
fy (eo = probability density function of D 


femene cumulative distribution function A. is the same for 


all members, 


ve ae | is FG]. f(x) dx (eqn 3.56) 


- © 
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where F, Ox) is the cumulative distribution function of A 
£5 (2x) Poetics puODaDIlityvedensity function of D 
As an approximation, an upper bound on Pe can be expressed 


as [Ref. 24]: 


N NN 
P.=1- [1 - 2 eee) at ree (XE. (x)..-..| £.(x)dx 
: } i=l 4 eee ry ? 
_ (ecinime sees 7) 
N N 
Ee: F. (x)f,(x)dx = = P 
| Sees izl 0 f 
rinen 
ae a aay < : P (eqn 3.38) 
: ti ls fexact i=] fy 
where Pee Zam noODaDiltty Of failure of the i-th member. 


7 
If the critical damage index and the design damage index 


memene Same for all members, then equations 3.36 and 3.3/7 
meance to: 


CO 


Pe cx [ F,(x)fp(x)dx = N Pe. (eqn 3.39) 


- © 


and 


P < P < N P (eqn 3.40) 
fi — fexact ~ 1 


where Pie = the probability of failure of any one member. 
1 
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IV. COMPUTATIONAL PROCEDURES 


Memevaluate the probability of failure for a given struc- 
fMimal component, @ constant amplitude sinusoidal S-N curve 
Meoation 3.13), statistical parameters of A, and either 
the spectral density function or a representative time history 
of stress response must be known. The overall computational 
Mmgcedire is shown in Figure 4.1, depicting the step-by-step 
procedure to calculate the probability of fatigue failure. 

In the following analysis, it is assumed that the spectral 


density response curve is available. 


A, DETERMINE CRITICAL DAMAGE INDEX SHAPE PARAMETERS 
For any Weibull distribution, the cumulative distribution 


@ame be defined as: 


n 
Box) = - - exp | = | (eqn 4.1) 


where & = Weibull scale parameter 


Il 


n Weibull shape parameter 

Memprobability plotting paper exists for direct plotting of 
distribution data which contains a shape parameter. However, 
for the Weibull distribution, the shape and scale parameters 
may be transformed to equivalent scale and location parameters, 


jms Weibull probability plots can be constructed [Ref. 16]. 


From equation 4.1, transformation yields: 
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n 
oo = argo || © ] (eqn 2) 
and then, 


Inin ET)! Selene x) aan & ) (eqn 4.3) 


Therefore, a plot of any Weibull variate Inin[1/{1 - F(X))] 
plots as a straight line against the natural logarithms of 
the observations. 
Rewriting equation 4.5 as: 
w= a+ bz (eqn 4.4) 


where W = Inln[1/[1 - F(x) ] 


ge In{x) 
by aa 
oe = TE CS ) 
then, 
n = b (eqn 4.5) 
and : 
5 = exp [ - 7 | (eqn 4.6) 


Where the values a and b are obtained as the y intercept and 
Slope of the plotted line. Within this study, b of equation 
fete is redesignated ¢ with appropriate subscripts (l.e., Cy 


or Cy) > and a of equation 4.4 is redesignated b with appro- 


priate subscripts (i.e., b, or Dy) . 
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INPUT DATA 


Constant Amplitude 
S-N Curve 


Specera lw Wensity. Funct ion 
G(w) of Stress Response 


Obtain Fatigue Data 









Statistically Reduce Datal——— >} Statistical Parameters 
Oe: 

Generate Probability 
Meer ibution Plot: 


Number of Points in 
Random Process Simulation 


Ultimate Strength 


Coefficient of Variation 
of Design Damage Index D 


Avs percent failure 


Simulate X(t) 


Evatuate No» Mora 


Transform X(t) to Peak and Trough Sequences of 
Stress Histories: Apply Rainflow Cycle 
Counting Method with Correction for Mean Stress 











Calculate A and o, and D and T, 


Determine Statistical Parameters of D: bp, cp 
Determine Central Safety Factor: n 


Evaluate Probability of Failure of Component 


Weibull Model Log-Normal Model 


ExXtmapolate System Pe from Component Pe 


Figure 4.] COMPUTATIONAL PROCEDURE FLOWCHART 
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PeeeeOUMULATE STRESS TIME HISTORY X(T) 
For a given spectral density function G(w), a sample of 


X(t) is simulated using the following expression [Ref. 25]: 


J 
mee) = 5 G(w-) Ow: cos(w.t + @.) (eqn 4.7) 
: 1 1 iL 1 
1=1 
and 
1 
7 = z Ws 5 Oc Ws < Wy (eqn 4.8) 
1=1 
where p. is a random phase angle, uniformly distributed in 
Ene Interval oO, 2m). 
Ow; — ene queneyerite Va | 
Ww, = frequency at the center of dw; 
me cffective frequency range in G(w) 
Mats Simulation of stress history is for a stress history 
with zero mean stress. Parametric studies of the effect of 


mean stress can be accomplished by using the relationship: 


X(t; ) = X(t) + THM (eqn 4.9) 


where THM is the mean stress amplitude to be studied. 


fe CALCULATE No» Mm, AND a 


When the spectral density function is known, n m 


Or Or 


meee Can be Calculated as follows: 
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l 2 w? G(w) dw 
ia = a= ae (eqn AO } 
a a fr G(w) dw 
and 
eS 4 
 . 1 [-s w’ Glw)dw 
O i S Meec curds (eqn 4.11) 


mem & can be calculated using equation 3.8. 


ieeeecALCULATE STRESS RANGES 

iranstorm X(t) to peak and trough time history and apply 
@memraintlow counting method to determine the stress ranges. 
In the developed source code, subroutine RNFLW counts stress 
ranges from troughs to peaks, and subroutine RNDRP counts 
stress ranges peaks to troughs. The mean and standard devi- 
fon of each stress range is then calculated by subroutine 
STAT. These statistical parameters are then averaged to 
obtain the overall mean and standard deviation of the simu- 


lated stress ranges. 


Eee fol IMATE Elie | 
ime expected value of s° with correction for the effects 
of mean stress can be estimated using one of the mean stress 


models previously discussed. For this study, the Goodman 
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B 


model was chosen. The expected value of S$ SOLRected 
for mean stress can be expressed as: 
B 
aye a 4.12 
Ip [S ] 7 i¢ ae Jo. (eqn sh ) 
i=l l ab 
-_=—— 
u 


where K = number of stress cycles counted by the rainflow 
cycle counting method 


oO. =mean stress of the i-th cycle 


F. CALCULATE A AND o, , AND D AND o 


A D 

From the Weibull parameters b, and c, obtained from the 
curve-fit procedured previously discussed, the mean, standard 
deviation, and the coefficient of variation for A can be deter- 
mined via equations 2.4, 2.5, and 2.6. As a result of the 
Tainflow cycle counting of the simulated stress response time 
history for the given mean and standard deviation of the design 


fatigue life, the mean and the standard deviation of the design 


damage index D can be calculated by [Ref. 23]: 


B = 
i 2 IS" ] Np (eqn 4.13) 
(on)? 
and 
E [s*] 
Ty = On (eqn 4.14) 
g D 
(o¢*) 
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tt 


where Ny mean design number of stress cycles 


"Ny Standard deviation of design number of cycles 


G. DETERMINE DESIGN DAMAGE INDEX SHAPE PARAMETERS 

The mean and standard deviation of the design damage index 
Previously calculated are fed to the subroutine WEIBR. This 
subroutine uses an iterative procedure on equation 2.6 with 
initial guess values of the shape parameters to generate a 
Meemticient of variation. This generated CV is compared to 
the actual CV determined from CV = op/D. The initial guess 
Values are then modified and the iteration repeated until the 
memeraced CV is within one percent of the actual CV. Once 
this criterion is met, the calculated shape parameters are 


Geeurned to the main program. 


femeecenLCULATE PROBABILITY OF FATIGUE FAILURE 
ieee Weibull Model 

Seno scoduartoneo. ZO, the probability of failure is cal- 
Sieaced Vila numerical integration. In that integration to 
Miminity 1S not possible, probabilistic theory of distribution 
Mmes 15 Utilized to establish the upper limit of integration. 
Beem probabilistic theory, it is accepted that 99.7% of the 
area under a distribution curve is contained within the range 
of the mean plus and minus three standard deviations of the mean. 

The numerical integration routine used, DQSF, uses the 


mean of either the design damage index, or the critical damage 


48 





index as the upper limit of integration for the first iteration. 
Miemeeper limit of integration is then incremented with one 
standard deviation of the appropriate damage index and the 
integral re-evaluated. This iterative procedure is continued 
fmeiecither the predicted probability of failure of the current 
iteration is increased less than one percent over the previous 
iteration, or ten standard deviations have been added. This 
procedure allows maximum penetration into the distributions 
fee region, and usually accounts for all but 6E-07 percent of 
the area under the distribution curve. 
2. Log-Normal Model 

Using equation 3.32 and the subroutine MDNOR, the prob- 

ability of survival is calculated. The probability of failure 


is then calculated as one minus the probability of survival. 


49 





V. NUMERICAL EXAMPLES 


A. INTRODUCTION 
Meorder to test the developed source code, a sensitivity 
analysis of several variables affecting the predicted prob- 
meerty Of failure was conducted. In that this study is a 
@mect expansion of published material on deterministic analysis 
fmeatigue life prediction [Ref. 6], the control variables were 
chosen to be the same as in that previous study, with the ex- 
Bemeron Of the time history data block duration. The original 
Study centered on fatigue life predictions for light water 
reactor components subjected to low amplitude, low frequency 
Stress fields induced by flow excitation. 
meee oucn opectral Density of Stress Response 

To simulate this low frequency, low amplitude stress 
response, a power spectral density (PSD) response function was 
assumed. The original study tested both a single peak and a 
multipeak PSD, and determined that a single peak PSD used to 
Simulate a multipeak PSD at the dominant frequency with equiva- 
lent RMS stress yields a conservative fatigue life prediction, 
thus is safe for use as an approximation. 

In that a single peak PSD can be safely used, all test 
runs were made using the same single peak PSD. The PSD used is 
depicted in Figure 5.1. The area encompassed within this simula- 


iron yields a Opus Ce wkoin wWolen IssOne Of Live values 
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previously tested. This PSD has a calculated irregularity 
meaeean © Of 0.935, implying that the PSD is a quasi-narrow 
Pemoedistribution, and that the Rayleigh distribution can be 
used to solve equation 3.4 for the value of the design damage 


index. 
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2. Time History DaAtawo foe kK 


Having chosen the PSD, equation 4.5 was solved for this 
mee, creating the simulated stress time, history depicted in 
mepure 5.2. in the original study, a time history data block 
of one second was used. For this research, a two second data 


block was chosen for the second analytic study. For the first 


Sel 





Bmmay, the time history data block duration was used as a con- 


trol variable, with the data block time duration ranging from 


Omemeo twenty seconds. 


The fundamental premise of a time history data block 


is that it repeats itself ad-infinitum, and thus, can be used 


bemoamiulate lives in the infinite lite region of the S-N curves. 


symmetry of the data block about its time axis midpoint insures 


meena tion of the block. 


The two second block appears to be 


more nearly symmetric than data blocks of a shorter time span. 
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TecOlurrmlechissObSenvartton, data blocks up to twenty 
seconds in duration were tested. It was found that even time 
blocks are essentially symmetric about their time axis midpoint, 
while the odd time blocks are not as close to being symmetric, 
thus their repetition time is longer than the time interval 
represented by the data block. None of the twenty examined data 
blocks exhibited exact symmetry, but all of the even blocks were 
essentially symmetric, and since the two second block had the 
tiemtcot repetition cycle, it was chosen as the time history 
data block duration cycle for use in the second study. Although 
being further from true symmetry than the even data blocks, the 
odd time data blocks were retained in the first study because all 
of the data generated as a result of using those data blocks fell 
femen the acceptance specifications of "Chauvenet's criterion" 
een 20; pp. 65-67]. 

fee Fixed Input Parameters 

The remaining input variables that were to be permanently 

mexca were taken from Shin's paper [Ref. 6], with the exception 


of the Weibull parameters DA and c which were determined from 


A b 
the data in Buch's paper [Ref. 14], as previously discussed. 
A summary of the input variables that were fixed throughout all 


testing is provided as Table I (variable names are fully defined 


in Appendix A). 
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a Densliavity Variables 


Micrreiameiny daha oko etter the nunbDer Of points in 
the random process simulation, the coefficient of variation 
Gieene design damage index and the number of stress cycles 
were the sensitivity variables. Additionally, the beta values 
were sensitivity variables, but the same five values were used 
Mmeeoushout the study, so they are listed as fixed parametric 
Mmmenese in the first analytic study, the coefficient of vari- 
meron was tested at 0.2, 0.5, and 0.8, and the design cyclic 


life was fixed at 1.8x1019 


cycles. This study was to evaluate 
the sensitivity of each model to the design damage index coef- 
mreient of variation. 


The second study was of the effect of cyclic life on 


the models. For coefficients of variation of 0.2 and 0.8, 


24 





be cycles, while the time 


cyclic life was varied from 10° to 10 
history data block was fixed at two seconds (i.e., DT = 0.002 
seconds, and IN = 1000). Plots of probability of failure vs 
cycles applied were generated for a single component, typically 
ranging from 10° ne iy Syeltcceecdaitionatiy=. plots of Lfirtty 
member structures were generated, with lives ranging from 10° 


jdt 


to 10 cycles, depending on the coefficient of variation 


Studied [Ref. 24]. 


Pemex eDUCTION OF BUCH'S DATA 

In order to test the prediction models, critical damage 
index parameters (1.e., mean, standard deviation, coefficient 
of variation and Weibull distribution shape factors) were 
required. Buch [Ref. 14] consolidated test data of several 
investigators, with the most comprehensive data available being 
for aluminum alloy 2024-T3 sheet specimens. Thirty-seven sets 
of these damage index values were reduced and plotted utilizing 
the methods presented in Chapter IV, section A, and shown in 
Figure 2.3. The Weibull parameters b and c were then obtained 
Via equations 4.5 and 4.6. 

It is recognized that the S-N curves for aluminum alloys 
Moe not have a classical "fatigue endurance limit" as do S-N 
Curves for some steels, thus the infinite life region is poorly 
defined. However, the applied PSD with a Tams Oe OP KS 1 aks 
below the approximate fatigue limit range for reversed bending 


Seealuminum alloys. The fatigue limit range is listed as $8.0 


a8 





Memmmoe oO Ksi [Reft. 27; p. 5-9], thus the "infinite life" assump- 
tion was deemed valid. Additionally, a review of the simulated 
relative stress response plot (Figure 5.2) shows that the maxi- 
mum stress level incurred is about 6.0 ksi, safely below the 
lower limit of the approximate fatigue limit range. 

Meenin Buch's paper, a limited amount of data on steels 
was presented. Although the data points were too few in number 
to obtain a representative data base for probabilistic reduc- 
tion, the discussed trends were the same as those for the alumi- 
Meeailoys. therefore, the specific distribution obtained from 
the aluminum data is representative of the statistical scatter 
for most metals, and was deemed acceptable for use as an illus- 
trative example. This observation was repeated in a review 
Oeeyersching’s work [Refs. 12,13] and [Refs. 28,29]. The trends 
of all materials studied, except a smooth specimen composite, 
Were virtually identical. 

The mean of the critical damage index for all of the mate- 
Yials was within the range of 0.836 to 1.72, excluding the com- 
Meoites, With coefficients of variation of approximately 0.65. 
This narrow band led Wirsching to use an assumed mean of 1.00 
pmomcy Of 0.65 for all probabilistic studies, regardless of 
the material involved. In that the precedent has been set for 
assuming a "typical" statistical distribution is valid for most 
materials, using fatigue strength exponents associated with 
Steels anda distribution formulated from failure of aluminum 


alloy samples was accepted as a reasonable combination of 
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parameters. This combination is possible only because of the 
wide range of statistical scatter in the available experimental 
data, and this approximation is presumed to be invalid for data 


sets having small coefficients of variation. 


eeeeeNoll TIVITY TO COEFFICIENTS OF VARIATION 

the coefficients of variation of the design damage index 
Studied were 0.2, 0.5 and 0.8. Gnas -22 ts representative of 
a distribution with relatively little statistical scatter, 
while CV, = 0.5 yields a Standard deviation of fifty percent 
of the mean value and represents a fairly large statistical 
eteccr. However, statistical scatter up to seventy percent 
is common, with scatter up to ninety-eight percent in rare 
cases, for fatigue studies [Ref. 28]. Thus, CV = 0.2 is 
deemed to be conservative, CV 5 = 0.5 is midrange, and CVy = 0.8 
1s representative of most "worst case'' scatter situations. 

The cyclic life values reported ranged between 1.5 and 
Me xX oe cycles [Ref. 23], for central safety factors between 
one and three. Actual values studied ranged from 10° 1) Lee 
maelmes, With central safety factors ranging from near zero to 
over one thousand. Common engineering safety factors are in 
mie range of one to three, frequently as high as five and occa- 
Sionally as high as ten, so the range of central safety factors 
from one to three received the most scrutiny, with some study 
of the range one to ten. Generated data outside this range 


was not plotted. The study used a control life of 1.8 x ge" 
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cycles, which corresponds to ten years continuous operation of 
a component subjected to the PSD specified. This value was 
mammnea trom the deterministic fatigue life prediction deter- 
mined by the deterministic portion of the source code. 

For each coefficient of variation studied, the time history 
data block was varied from one to twenty seconds by changing 
the number of points in the random process simulation. In 
m@enrements of 500 points, the point process simulation was 
varied from 500 to 10000 points. The time increment was fixed 
at 0.002 seconds per point. By varying the time data block, 
the stress response function X(t) was changed. This, in turn, 
changed the number of simulated cycles counted by the rainflow 
counting method, thus E[s?} and D were varied. 

For each value of D, the corresponding central safety fac- 
meme equation 3.25) and the predicted probability of failure 
fOr C€ach probability distribution model were calculated. 
Graphical presentations of the probability of failure versus 
m@mes central safety factor for each distribution, separately 
eo together, were generated. 

As can be seen in Figure 5.3, the Weibull model is sensi- 
mye to the coefficient to variation up to a central safety 
meteor value of about 2.25. Beyond this value of n, all three 
CV's studied converge and predict the same probability of 
failure. Any further change in the time history data block 


duration has no affect on the separation. Although not plotted, 
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Seeaeruns for central safety factors between ten and twenty 
were conducted, with no separation in predicted values. 

Maaco +eresents the results of the log-normal distri- 
Piition. The log-normal model is sensitive to the coefficient 
Seevertation throughout the entire range of central safety 
Pieters studied. Below a central safety factor value of about 
2, CV, = 0.2 is the most conservative and CV 5 =—Omoews. tae 
least conservative. At a central safety factor value of about 
Eyeeene predicted values cross over, with CVp = 0.8 being the 
most conservative, and CV, = 0.2 Deing the least conservative. 
memeall values of n, CVp = 0.5 tracks between the limits estab- 
[rienced by the other two coefficients of variation. 

When the Weibull and log-normal models are compared, each 
predicts approximately the same probability of failure for 
M@emitral safety factors near one. This compatibility of models 
was expected, and was one of the acceptance criteria for the 
feemeee Code developed. in the central regions of fatigue data 
distributions, the log-normal and Weibull distributions are 
Peren indistinguishable [Ref. 30}. In this specific case, 
ieee 1.0 means D = A , and A was fixed at 0.9, the mean of the 
distribution being studied. 

padi tionally, there 1s almost always a significant differ- 
ence in the predicted values within the tail regions of the 
distribution, with as much as a one order Magnitude Gditrerence 
when the actual value is twenty percent or more below the mean 


Pemue Of the index [Ref. 30]. Figures 5.5, 5.6, and 5.7 
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demonstrate this trend is exhibited by the developed code pre- 
dictions for each coefficient of variation studied. The 
Seeameted probability of failure is nearly equal for the two 
Meaers for a central safety factor of one, quickly diverge, 
and then settle out with a near constant separation. Figure 
5.8 is included to demonstrate that this separation exists 
well into the tail region, but decreases as the elapsed life 
becomes smaller in magnitude. 

Mmeeiat the variable component of the central safety factor 
is the cumulative design damage index D, and D is a function of 
eiapsed cycles, the higher the value of n, the earlier in 
Bamuoue life the model is analyzing. Therefore, the prediction 
models meet the expected separation of predicted values in the 
memieGegion criteria established by the models predecessors. 

Based on the coefficient of variation sensitivity analysis, 
the developed source code does what it was expected to do. 

That is, it generates essentially equivalent predicted prob- 
abilities of failure for the Weibull and log-normal model in 

the central region of the distribution, and demonstrates the 
expected divergence of predicted values as the tail regions 

are entered. In all reference material reviewed, the log- 
normal model is considerably less conservative than the Weibull 
Megel in the tail regions of any distribution, whether the orig- 


meciedata tit the Weibull or log-normal distribution better. 
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WEIBULL vs LOG-NORMAL 
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WETBULL vs LOG-NORMAL 
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fees NotTIVITY TO APPLIED CYCLES: SINGLE MEMBER 

Having examined the models sensitivity to the coefficient 
Pevariation of the design damage index, testing proceeded, using 
applied design cycles as the sensitivity variable. Based on the 
results of the first test phase, the time history data block 
was fixed at two seconds (i.e., 0.002 second time interval, with 
1000 points in the random process simulation). Additionally, 
the design damage index coefficient of variation was only tested 
at CV) eeeeze and 0.5, In that the Weibull model showed conver- 
aomecmor predicted probabilities of failure for high central 
safety factors (i.e., low magnitudes of applied cycles), with 
CV, = 0.5 tracking between CV, =o comand. in tne Sensitive 
range, and although the log-normal model was sensitive to CV, 
am all ranges, CV) = 0.5 tracked consistently between the 
bounds of the other two test values, even after cross over 
occurred, the second phase testing was conducted without CV), = 
Was a test value. 

femiimetme first study, design cyclic life was 1.8 ioe 
mates. Beta values of 12, 12.5, 13, 13.5 and 14 were tested 
for each design damage index coefficient of variation and applied 
feess cycle value tested. For single member testing, applied 
cyclic life was varied from i0° to ie CG lose 

iesueero.9 1S a presentation of predicted probability of 
failure versus number of cycles for the Weibull model, for 
CV, = 0.2. This plot shows that for B = 12 (the structure with 


the lowest resistance to fatigue failure), the probability of 
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mBaelure is about one hundred percent at 1.8 nO eye ites. ehele! 


as the beta value is increased (i.e., higher resistance to 
fatigue failure), the components cyclic life is increased. 
As beta is increased linearly, the design fatigue life in- 
creases exponentially. This is amplified in Figure 5.10. 

When the data of Figure 5.9 is plotted on a log-log scale, 
as in Figure 5.10, when 6 = 12, the probability of failure 


at 1.8 101° 


cycles is one hundred percent, while for ® = 14, 
the probability of failure is only about two percent. 

moure 5.11 is a presentation of the predicted probability 
of failure versus number of cycles for the log-normal model, 
fOr CV, = 0.2. This plot exhibits the same shape as the 
Weibull plot, but in the low probability of failure region, 
the log-normal model predicts an approximate one order of 
magnitude higher cyclic life before component failure begins, 
memc slightly lower cyclic life to one hundred percent fail- 
ure. These differences are amplified in Figure 5.12, which 
depicts a comparison of the two models. As in the coefficient 
of variation study, the models are virtually indistinguishable 
in the central region, and are about one order of magnitude 
Meare in the left tail region of the distribution. As pre- 
viously discussed, these observations were expected, and felt 
Membe a necessary occurrence in order to establish the source 
code methodology as valid. 


Figure 5.13 depicts the log-normal model on a log-log plot, 


emphazing the left tail region performance. Figure 5.14 compares 
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the performance of the Weibull and log-normal distribution 
Mogers in the left tail region. Again, the one order of mag- 
nitude separation is exhibited in the left tail region, and 
Mear identical predictions in the central region of the 

aestr bution. 

imeene log-normal plot, the structure with the lowest re- 
mreimee tO fatigue failure (i.e., 8 = 12) reaches one hundred 
percent predicted failure slightly befor 1.8 ie cycles. This 
is due to the poor resolution of any log-normal probability 
model deep in the tail region of a data distribution. This ob- 
servation is repeated for all beta values. That is, complete 
failure is first predicted by the log-normal model, and onset 
of failure is first predicted by the Weibull model. 

To verify that the established trends were applicable for 
Mmewer values of the coefficient of variation of the design 
damage index, the analysis was repeated for CV) = eo ouGes 
5.15 and 5.16 show the Weibull and log-normal models performance 
iM~mime same manor as Figures 5.9 and 5.11 did for CV, = Oe 
Beoure 5.17 is a comparison of the two models for CV, = 0.8. 

In all three plots, the shape of the curves and trends of the 
plots follow the same pattern established by the CV, = Gee 
BELOtS . 

Pecume os. 18 1S a comparison of the coefficient of variation 
Meiisitivity as a function of applied cycles for the Weibull 
model. As was the case in the CV, sensitivity study, for low 


cyclic lives, the model is insensitive to CV, with negligible 
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Taeweicron in predicted failure probabilities. As cyclic 

life nears design life, there 1s a separation in predicted 
values, with the CVy = 0.2 being the more conservative value. 
From the source code generated data, the CVn = Ue plots do 
ewenemally reach a predicted probability of failure of one 
Mgm@e@nred percent, but so far into cyclic life, that it does not 
Show on this plot. 

Figure 5.19 is a comparison of the coefficient of variation 
wememcrvyity for the log-normal model. Again, the predicted 
trends are sensitive to the CV 5 throughout the range of study, 
tmaiea noticable cross over point. Early in cyclic life, 

a 0.8 is the more conservative value, and near design 


Byte life, CV, = 0.2 is the more conservative value. 


Meee NSETIVITY TO APPLIED CYCLES: MULTIMEMBER 

Fatigue life of a weakest link structure with fifty members 
was evaluated. Again, beta values of 12, 12.5, 13, 13.5 and 
14, in conjunction with design damage index coefficients of 
Warlation of 0.2 and 0.8, with a fixed time history data block 
duration of two seconds, were tested. All members were assumed 
to follow the same probability model for the damage index, and 
all members were specified to have the same probability of 
failure for a given loading spectrum and elapsed stress cycles. 

The upper bound for failure of the structure was determined 


and plotted in terms of total number of applied cycles, using 
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the procedure and equations outlined in Chapter ii esecetion G, 
For a design damage index coefficient of variation of 0.2, 
Figure 5.20 presents the results of the Weibull model, and 
Figure 5.21 presents the results of the log-normal model. 
Figure 5.22 compares the two models, with the Weibull model 
being significantly more conservative than the log-normal model. 

As Figure 5.22 demonstrates, very few additional cycles can 
be placed on a log-normal modeled system before predicted fail- 
ure transits from near zero to one hundred percent. The Weibull 
model, on the other hand, does have an obvious transition range. 
The presence of this transition range 1s very beneficial if 
Maintenance and inspection procedures are incorporated into 
Paetouc life prediction modeling. 

mmeoures 2.25 through 5.25 present the same information for 


CV, = 0.8. Figure 5.26 compares the Weibull model for the two 


D 
CV's. As expected, there is a negligible separation in the 


D 

curves for each beta value, due to insensitivity to CV 5 at low 
Meme lite, Figure 5727 compares the log-normal model for 
igae two CV's Studied. This plot does show a Separation in 
predicted value for each beta, for each CVy> When 1S due Eo 
the log-normal model being sensitive to the CV, throughout 
the range of analysis. 

Figure 5.28 is a comparison of the Weibull and log-normal 
models for a combination of single and multimember predicted 


results. In that for the entire range of beta's studied, the 


pots all have the same shape, a single beta value is presented. 
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Figure 5.20 WEIBULL MULTIMEMBER LIFE CYCLE SENSITIVITY CVp = 0.2 
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Figure 5.21 LOG-NORMAL MULTIMEMBER LIFE CYCLE SENSITIVITY 
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WEIBULL vs LOG-NORMAL 
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Figure 5.22 WEIBULL VS LOG-NORMAL MULTIMEMBER COMPARISON 
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WEIBULL vs LOG-NORMAL 
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Figure 5.25 WEIBULL VS LOG-NORMAL MULTIMEMBER COMPARISON 
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LOG-NORMAL MODEL 
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Figure 5.27 LOG-NORMAL MULTIMEMBER CVp = 0.2 VS CVy = 0.8 
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Mmenis case, beta equals 153 was chosen to represent all beta 
Memes. A review of Figure 5.28 shows a Significant separa- 
tion in single member and multimember predicted life for the 
Weibull model. The observed separation in predicted life is 
somewhat more than a complete order of magnitude in cyclic 
life. For the log-normal model, the separation in predicted 
life is only about one-half of an order of magnitude. For a 
system in which the probability of failure is the same for all 
components of the system, at least a one order of magnitude 
difference in a single components life and the system life 
would be expected. Therefore, the Weibull mcdel is more rep- 
Pesentative of a systems extrapolated fatigue life. This is 
due primarily, to the response of the Weibull model in the 
meee tail region of the distribution. 

This adaptability of the Weibull model to system life 
extrapolation was expected, in that the Weibull distribution 
has been extensively studied in its application to system lives 
[Ref. 31]. Other probabilistic distributions have been adapted 
to the weakest link analogy, but the Weibull distribution is 
Stili considered the best for this purpose. Many investigators 
argue that the Weibull model is too conservative in the left 
meals region of a fatigue data distribution (1.e., Wirsching, 
et. al.), but proponents of either distribution model do concede 
that the Weibull model is best for system life extrapolation 


met. 30]. 
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VI. SUMMARY AND CONCLUSIONS 


A. SUMMARY 

Peorobability based approach to evaluate the high-cycle 
fatigue life of single component and multiple member, weakest 
link structures, in random stress environments, was described. 
iieceapproach utilizes the Palmgren-Miner linear damage rule, 
with compensation for mean stress effects and the statistical 
variability of the critical damage index. The cumulative design 
damage index D was assumed to be a random variable, following 
the same probability distribution as the critical damage index 
feeeecor the particular probability foam ieloOnmebelms analyZed 
(i.e., Weibull or log-normal). Cyclic stress history was 
assumed to be representable by a finite, repeating, time history 
data block. 

Srauistical reduction and best fit analysis of selected 
existing experimental data for A was accomplished. The subse- 
meen application of ordér statistics and probability paper 
ouemeeor the data produced a best fit Weibull distribution, 
Teed mean of 0.90 and a coefficient of variation of 0.6/7, 
Vice the proposed model [Ref. 13], of an assumed log-normal 
rst Gipution with a mean of 1.00 and a coefficient of varia- 
tion of 0.65. The calculated Weibull distribution shape para- 


meters by and c, were Oey Aiaiceele o2e) respective ly. 





Hog-normal models can be formulated as a function of the 
mean and coefficient of variation of the data distribution, 
ferme are the same for either log-normal or Weibull distribu- 
meee. iheretore, the developed Weibull model was tested by 
comparing it to the values of fatigue life obtained by assunm- 
memene data did fit a log-normal distribution, and solving 
the closed form mathematical representation of the log-normal 
ret? tDuUtion. 

Bensitivity studies of the predicted probability of fail- 
memes dependent on resistance to fatigue failure ({i.e., nega- 
tive reciprocal of the fatigue strength exponent), coefficient 
@mevariation of the cumulative design damage index, central 
mimety factor, and applied stress cycles were conducted. For 
both probability distribution models, as beta increases linearly, 
Emme life increases exponentially. The log-normal model is 
Sensitive to the coefficient of variation of the design damage 
index throughout the entire range of analysis (central safety 
factor or applied stress cycles), but the more conservative 


feeebowo. CV 15 dependent on elapsed cyclic life, with a pro- 


D 
mommecd €ross-over. The Weibull model is insensitive to CV 
mimearily Cyclic life (i.e., elapsed life less than about eighty 
percent of design life). During that portion of life in which 


the Weibull model is sensitive to CV the lower the value of 


le? 
CV ps the more conservative the predicted fatigue life (i.e., 
Mmieemicher the predicted probability of failure for a given 


value of the sensitivity variable). 
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for all selected methods of graphical presentation of 
results, the predicted probability of failure for the two 
models yield comparative observations that agree with the 
published behavior [Refs. 10, 13, 30, and 32] of the two dis- 
tributions. That is, in the middle region of the data distri- 
bution, the two models are virtually indistinguishable. In 
the left tail region of the data distribution, the Weibull 
model is significantly more conservative than the log-normal 
model, with up to a one order of magnitude difference in pre- 
dicted fatigue life. 

When the weakest link analogy was used to extrapolate an 
upper bound of a multiple member structures fatigue life, the 
Weibull model was considerably more conservative in predicted 
fatigue life. The separation in the two models ranged from 
about one order of magnitude at end of life to two orders of 


magnitude in the left tail region. 


Pee CONCLUSIONS 

Extensive investigation of, and efforts to develop prob- 
foietstic fatigue life prediction models has occurred. Early 
models were developed using cyclic life as the basis of 
analysis, due to the availability of cyclic test data. In 
low cycle fatigue regimes, test data can be produced in a rela- 
tively short time span. In the high cycle fatigue regime, ex- 
tended periods of testing are required, but test data can be 


obtained. In the ultra-high cycle fatigue regime, it can 
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take years to cycle a single sample the required number of 
times to represent the components life. To get a representa- 
tive data base, testing costs can exceed component life-cycle 
costs, thus prediction models that can use a much more limited 
data base are desirable. 

From the above requirement, fatigue life prediction models 
using the Palmgren-Miner linear damage law and known statistical 
distributions of fatigue data are being developed. A FORTRAN 
computer code which utilizes this method was developed and 
tested, using both the Weibull and log-normal distribution models. 
Testing of the code demonstrates that the generated data complies 
with the published trends of prediction models formulated on 
cyclic life for these two distributions. 

New structural designs entail significant modeling and/or 
testing. Components designed to operate in the ultra-high 


meee iatigue regime can entail cyclic testing up to ia 


Gyveles 
fOr a single sample, and then repeating the test on "identical" 
Samples a sufficient number of times to get a representative 
data base. This model can be used to selectively reject mate- 
Yial compositions which appear to be feasible, but probabilistic 
analysis imply will fail earlier than design requirements allow. 
Having used the model as a first cut analysis of potential 
design components, testing expense can be reserved for viable 
compositions, not spent on compositions that might work. 


An extensive literature search has demonstrated that there 


are as many proponents for using the Weibull model as for using 
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the log-normal model. The log-normal model is easy to formu- 
late, and has extensive library routines providing closed 

form solution techniques. Its simplicity of implementation 
makes it a very popular model to use, and more importantly, 
many existing sets of fatigue data fit the log-normal distri- 
bution. The Weibull model is not as easily implemented. The 
Weibull shape parameters must be calculated, and then the mathe- 
matical formulation must be numerically integrated. The major 
advantage of the Weibull model is its adaptability to multiple 
member structures fatigue life extrapolation from single com- 
ponent predicted fatigue life values. 

The intention of this study was not to prove the Weibull 
model superior to the log-normal model, or vice versa. The 
wemieered experimental data fit the Weibull distribution better 
moment fit the log-normal distribution, so it was known that 
the Weibull model would be superior in this case. Many .data 
Sememecan be found in which the log-normal distribution fits 
Memeer, thus the log-normal distribution would be superior. 
itespurpose for using both distributions was to have a conm- 
parison method for analyzing the developed Weibull models 
predictions. The comparative testing produced observations 
consistent with all comparisons located in the literature, 
thus implying the formulation of the models was correctly 
eecomp lished. 

Paccamom the mmterrelationships of the two models, and 


the compliance with published behavior of existing fatigue 


oF 





meecmprediction models, it is felt that the developed source 


code is an acceptable and viable fatigue life prediction tool. 
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VII. RECOMMENDATIONS FOR FUTURE WORK 


A. GENERAL RECOMMENDATIONS 

In that the developed source code is new and was not drafted 
by a programmer, some modifications to the code which enhance 
Eeemerency of computer cost factors are possible. A significant 
cost factor is the numerical integration routine used. There 
ieee more efficient and/or more accurate library routines 
available for numerical integration. Additionally, as submitted 
in Appendix B, the source code contains numerous write and asso- 
ciated format statements which were beneficial in code imple- 
mentation, but may be deemed superfluous information in general 
ce. 

The code was tested using only one set of fatigue data, and 
more importantly, associated control parameters were not consist- 
Siewith the material properties of the source data. It is recom- 
mended that additional data be obtained, and a rigid adherence 
of consistent material property parameters be applied to the 
feweadta to confirm predicted probabilistic fatigue life complies 
with the data sets life. Also, it may be a beneficial exercise 
memsoelect a data set that fits the log-normal distribution better 
than it does the Weibull distribution, to analyze it and compare 
the results of that study with this study to evaluate model 


performance based on statistical fit of the data. 
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Most probabilistic fatigue life prediction models utilize 
the two parameter Weibull or log-normal distribution. In an 
effort to enhance the precision of cyclic life models, recent 
Studies into the use of the three parameter distribution has 
been conducted [Refs. 10,32]. In these studies, the third 
parameter is called the least life parameter, i.e., minimum 
cycles to failure. In terms of the critical damage index, a 
least life parameter would be the value of the cumulative 
damage index at which first failure occurs (assumed to be 
zero in the two parameter model). For initial implementation, 
the least life parameter could be assigned the value of A in 
the selected data set with the lowest magnitude. The selected 
data set would have provided a least life A of about 0.01, but 
from a review of the literature [Ref. 12], the least life value 
1s more generally around 0.18 for log-normal distributions. 

Using a least life parameter may reduce the left tail 
region separation of the two models. Additionally, incorpora- 
tion of a least life parameter could be beneficial, not only 
in enhancing accuracy of predicted life compared to actual 
life from test data, but may ease the implementation of main- 


tenance and inspection procedure effects analysis on the model. 


B. MAINTENANCE AND INSPECTION 
Reliability analysis of structures subjected to scheduled 
maintenance and inspection procedures has been investigated, 


primarily by engineers in the aircraft industry. Efforts to 
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incorporate the probabilistic analysis of crack detection, and 
subsequent repair or replacement of the faulty component, 
effects on structural fatigue life are being pursued [Refs. 
33,34]. When a crack is initiated, the components fatigue 
strength progressively decreases as the crack propagates, in- 
creasing fatigue failure rate with time. If the faulty com- 
ponent is detected by inspection and repaired or replaced, the 
Static and fatigue strength of the structure are renewed. This 
regeneration of component fatigue strength enhances the struc- 
ture's fatigue resistance and life. 

One possible way to account for the effects of inspection 
and maintenance procedures would be to adjust the least life 
parameter upward, commensurate with the reliability of the 
maintenance and inspection procedure. In that detection, thus 
memgearr of a crack is a probabilistic function of flaw size and 
resolution of the inspection procedure used (i.e., visual in- 
spection, liquid penetrant, ultrasonic, magnetic particle, 
Magnetic field perturbation, or radio-graphic), some cracks are 
not detected, so fatigue life does not restart after inspection. 
However most cracks that would cause early failure, represented 
by the left tail region, would be detected, and eliminated [Ref. 
S44). The predicted fatigue life can then be recalculated, re- 
jecting those sets of A deep in the left tail. Or, if the 
least life parameter were used, the tail region could be re- 
tained, with the least life parameter magnitude adjusted upward. 


ise allows use of the entire statistical distribution, and 
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accounts for the increased probability of survival. Rejecting 
data elements in the tail region requires recalculation of the 
Meremeana coefficient of variation of the “new” distribution. 
If the Weibull distribution were used, recalculation of the 
shape parameters is also required. 

Bperuectures other than aircraft undergo scheduled maintenance 
and inspection procedures (i.e., offshore structures, bridges, 
pimeidings etc.). Thus, incorporation of accountability for 
these procedures would enhance the usefulness of any prob- 
mommbistic fatigue life prediction model. Therefore, investiga- 
mieieincto incorporation of maintenance and inspection procedures 
accountability in the developed fatigue life prediction model 
is the primary recommendation for future work. It is believed 
Maatemodifying the source code to use the three parameter vice 


the two parameter distributions would ease this implementation. 


Seu LTIAXIAL FATIGUE 

A literature search was unproductive in locating material 
Bevering probabilistic fatigue life prediction for components 
Pwipmected to multiaxial fatigue. Extensive research is in 
progress to evaluate multiaxial fatigue in the deterministic 
sense, but no consensus yet exists as to one method of analysis. 
Numerous proposals for high cycle fatigue analysis have been 
generated. Garud [Ref. 35] presented an excellent state-of- 
the-art survey of multiaxial fatigue and discusses most of the 


proposed models. 
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Each proposed model appears to be limited to a specific 
ipeemeation, 1.e€., proportional cyclic loading in all axial 
@mrections, in-phase loading or out-of-phase loading, etc. 
Until an accepted correlation is developed, complex multi- 
axlal fatigue analysis of predicted fatigue life must be 
deferred well into the future. On the other hand, simple 
cases Of biaxial fatigue may be solvable now. In the high 
cycle fatigue regime, equivalent stress, determined by the 
Tresca or Von-Mises criterion can be used [Refs. 10,35,36]. 
However, recent studies indicate this procedure only applies 
for in-phase loading with fixed principle axes [Ref. 35]. 

For this limited case, investigation into using a repre- 
Sentative PSD with the rainflow counting method to generate 
a mean stress for each axis, then application of either the 
Tresca or Von-Mises criterion to generate an equivalent stress 
Can be accomplished. This equivalent stress would then be 
used in evaluating E[S®], with the remaining analysis pro- 
Beeamne aS Outlined in Figure 4.1. The predicted failure 
(1.e., elapsed cycles to fifty percent failure) could then 
be compared with the cycles to failure determined via the 
Classical deterministic approach of S-N curve analysis 
(outlined in reference 36). 

An additional study could be made by developing a PSD 
representative of the equivalent stress from the Tresca or 
Von-Mises criterion, then using this PSD with the rainflow 


cycle counting method as if the study were initially uniaxial. 
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iieu, comparison of application of equivalent stress before 
Sete Stress Simulation and application of equivalent stress 


meeet cyclic stress simulation could be conducted. 


D. SUMMARY 

The growth potential of the developed source code is un- 
limited. Some general improvements have been discussed, and 
two in-depth expansions have been proposed. Based on background 
material available, incorporation of maintenance and inspection 
procedure effects on prolonging fatigue life appears to the 


Subject area with the most promising, immediate, applicability. 
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Pies ND LA 


COMPUTER CODE USERS GUIDE 


ee NI RODUCTION 

The computer code ''FATIGUE FORTRAN" is designed to cal- 
culate the deterministic fatigue life, Weibull model prob- 
ability of failure and the log-normal model probability of 
memeore fOr a single structural component. Any one, or a 
Selected combination of the three models can be chosen. 
Peewtionally, calculations accounting for, or neglecting, 


mean stress effects are possible. 


BeeeeoOMPUTER CODE DESCRIPTION 

Representative time history data blocks, applied power 
Spectral density functions, material fatigue strength expo- 
hents, Weibull model shape parameters for the critical damage 
index, coefficient of variation for the design damage index 
and design cycle life to be analyzed are inserted as variables. 
The program then generates a random version spectral density 
function and converts the time history data block to relative 
stress amplitudes for the specified time sequence. The simu- 
lated stress history is then changed to a peak and trough rep- 
resentation for Rainflow cycle counting and subsequent calculation 
of the mean and standard deviation of the simulated cyclic 
Beness history. The deterministic fatigue life is then cal- 


culated for each beta value to be analyzed. 
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The program then calculates the Weibull shape parameters 
for the design damage index from the mean and standard devia- 
tion of the stress history simulations, as well as the mean, 
Seam@eand deviation and coefficient of variation of the cri- 
tical damage index from the critical damage index shape para- 
meters provided as input values. From these values, the 
Weibull model distribution and density functions are calculated 
and inserted into a numerical integration routine for deter- 
femacion of the probability of failure. The integration 
routine used is DQSF, a double precision routine using a com- 
feteedon Of Simpson's rule and Newton's three-eighths rule. 

The Weibull model formulation implies integration to in- 
jmommety, which is not realistic. To alleviate this problem, 
miemealling program is established so that the first iteration 
of the integral is evaluated from zero to the calculated mean 
of the damage index, (both the design and the critical damage 
index are used). Subsequent iterations are conducted using 
the damage index with one standard deviation of the damage 
igliex added for each iteration. From probabilistic theory, 
it is accepted that the mean plus three standard deviations 
Seeemeemean for a given distribution will include 99.7% of all 
Palements of that distribution. Rather than limit the integra- 
mitenmeto this level, the iteration is conducted until the dif- 
ference between successive iterations is less than one percent, 


Or ten standard deviations have been added. This typically 
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results in five standard deviations being added, with resul- 
tant inclusion of all but 6E-07 percent of the area under the 
Miseribution curve. 

The Weibull model is evaluated using both the critical 
damage index and the design damage index with their associated 
standard deviations as the upper limit of integration. The 
two methods yield essentially the same predicted probability 
Seeecariure. For parametric studies of the critical damage 
index, the user can chose to comment out the design damage 
index analysis section. 

The program then proceeds to the log-normal model, with 
Mmeemcirst step being the calculation of the log-normal shape 
factor. This shape factor is referred to as the "central 
safety factor,'' and is the ratio of the means of the two 
damage indexes. The parameters of the standardized argument 
for normal distributions are then evaluated and the standard- 
ized argument is then calculated and passed to a routine called 
MDNOR. This routine then returns the probability of survival 
Gometne specified standardized argument. The probability of 
fee 1s then calculated as one minus the probability of 


Survival. 


feeeecoOMPUTER CODE SUBROUTINES 

Three library routines are used, and are imbedded in the 
attached program listing. The MDNOR and GGUBS routines are 
Poom the IMSLDP library and are precompiled, so imbedding is 


not required. The DQSF routine is from the NONIMSL library 
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and 1S also precompiled, thus imbedding is not necessary. 
The remaining imbedded subroutines are non-library routines 
and are tailored to this particular program. Table II pro- 
vides a summary of subroutine names and purpose. 

The only subroutines with internal write statements are 
See@e@leand WEIBR. In both cases, for parametric studies, these 
write statements and formats can be commented out if desired. 
SPECT generates the relative spectral density parameters, and 
these values do not change for a given PSD and time history. 
WEIBR generates the Weibull shape factors for the design 
damage index. These values do change during parametric 
Studies if the coefficient of variation of the design damage 
index is one of the varied parameters. However, the calculated 
values are returned to the main program for use in the Weibull 


model section, thus writing these values is not necessary. 


TABLE Sit 


COMPUTER CODE SUBROUTINES 


RANSPC generates random spectral density function 
GGUBS basic uniform random number generator 

STAT evaluates mean and standard deviation 
RNFLW Rae low (eye lemecoune tiermcrough to peak 
RNDRP Paint lew Gy c¢ lemcountanowmpcak to trough 
SPECT computes standard deviation, irregularity 


factor, expected rate of zero crossings 
and expected rate of peaks 


WEIBR calculate Weibull parameters b and c from 
the mean and standard deviation 

DQSF numerical integration routine using Simpson's 
rule and Newton's three-eighths rule 

MDNOR Normal or Gaussian probability distribution 


fuincte Lom evaluator 
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D. COMPUTER CODE INPUT PARAMETERS 


A sample input file is presented as Table III for demon- 


peeamion Of format specified location and argument type, i.e., 


Character string, real or integer value. A description of 


each variable as encountered in the input file follows. Each 


parameter is referred to as it is named in the program. 


First Row 


1) 


Second 


it) 


2) 


3) 


2) 


>) 


6) 


7) 
8) 


me: 


Row 


NP: 


EN: 


NE: 


Da: 


DE: 


THM: 


Sus 


CiramadecteuectwinGmtttle sOL Specific test “case. 


Number of input frequencies of stress response 
power spectral density (PSD) function. 

Number of points in random time history 
Spunulllehecetey et: 

Number of BETA values (negative reciprocal of 
fatigue strength exponent) to be considered. 
Time interval in random time history simulation 
(sec). 

Average frequency interval to characterize the 
power spectral density function (Hz). 

Mean of time history relative stress response. 


Ultimate strength of component. 


IFLAG1: Integer flag value to specify type of analysis. 


=-Q0 Deterministic Fatigue Life analysis 
=] Probabilistic Fatigue analysis 


=7  oyeie}a 
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iiird Row 


1-5) 
Fourth Row 
i; 5) 
24,0) 
Fifth Row 
1) IWL: 
me) BCL: 
eB) CU: 
mm) BDEL: 
fo) mC DEL: 
Oo) MECVWI: 
On er CL: 


EX? ip DeiwsaGnay. values Of BETA to be analyzed. 


WF(NP): Frequencies of input PSD. 


SP(NE): Amplitude of PSD at specified frequency. 


Integer flag value to specify type of prob- 
abilistic analysis. 

=0 Weibull model 

=] log-normal model 

=2 both 

Lower limit of variability search for use in 
determining Weibull shape parameters for the 
design damage index. 

Dppouw iE wOne wat laud Llty Search - 

Weibull shape parameter for critical damage 
index (intercept of Weibull critical damage 
index: distri buc zon): 

Weibull shape parameter for critical damage 
index (slope of Weibull critical damage index 
Gis Ga lt home me 

Coefficient of variation of the design damage 
Tiga ee 


Design life cycles to be analyzed. 
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E. VARIABLE USAGE 

Persp. time history data block, and material (ultimate 
strength and fatigue strength exponents) are chosen and fixed. 
Then parametric studies are accomplished by varying the design 
damage index coefficient of variation and cyclic life. Depend- 
ing on the chosen coefficient of variation, the upper and lower 
limits of variability search may have to be changed. For 


0.1 < CVWM < 1.0, the resident values of CL and CU are adequate. 


fee COMPUTER CODE USAGE 

The program is tailored for usage on the IBM 3033, with the 
Bentix compilier. Unit five is used for input, and unit six is 
used for output. In addition, unit seven is used for creating 
a data file for subsequent graphical presentation of a time 
meoeory plot, and unit eight is used for creating a data file 
for subsequent graphical presentation of a simulated power spec- 
Mmermedensity plot. The data from unit seven can be placed 
meercecly into a plot routine, such as PLOTG, and submitted to 
mies versatec plotter. The data on the unit eight file must be 
used as input to a Fast Fourier Transform routine, which gener- 
mees tne data points for use in a plot routine. 

EOuweuarametric Studies, it 1s recommended that after the 
initial data run, the write statements and associated format 
statements for the unit seven and unit eight data be commented 
out. Additionally, a significant number of the unit six write 
Statements and associated format statements can be commented out, 


leaving only the required information for the particular study. 
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Program array dimension sizes dre listed on the first 
page of the program. Recommended values, as functions of 


the input variables are provided. 
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